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We show through rigorous calculations that dielectric microspheres can be organized by an incident 
electromagnetic plane wave into stable cluster configurations, which we call photonic molecules. The 
long-range optical binding force arises from multiple scattering between the spheres. A photonic 
molecule can exhibit a multiplicity of distinct geometries, including quasicrystal-like configurations, 
with exotic dynamics. Linear stability analysis and dynamical simulations show that the equilibrium 
configurations can correspond with either stable or a type of quasi-stable states exhibiting periodic 
particle motion in the presence of frictional dissipation. 
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In an optical field of non-uniform intensity, electromag- 
netic forces can move small particles towards regions of 
high intensity, and such optical gradient forces Q, 13, Q 
have been employed fruitfully to manipulate small parti- 
cles 0, 0] . Burns et. al. have demonstrated the exis- 
tence of another weaker optical binding force between 
two dielectric particles, induced by the multiple scat- 
tering. In contrast to van der Waals force, the optical 
force can be either attractive or repulsive. While Burns 
et. al. have already demonstrated the existence of such 
a force between two particles, and provided an under- 
standing within a simplified dipole picture, we here show 
through rigorous calculations that such kind of forces 
can organize a collection of particles into stable equilib- 
rium configurations that behave like "molecules" . The 
photonic molecules have scale length corresponds to the 
light wavelength, and the binding comes from multiple- 
scattering induced long range force that is strong enough 
to overcome van der Waals forces and thermal fluctua- 
tions. The same number of particles can be stabilised in 
a large variety of distinct shapes and bond lengths, and 
exhibit a multiplicity of static and drifting equilibrium 
configurations. What is rather amazing is that the clus- 
ters can have distinct shapes and structures, and have 
well-defined vibration frequencies, and yet the structural 
order is derived from an incident electromagnetic plane 
wave that has uniform intensity (in contrast to gradient 
force that this derived from a non- uniform field). Since 
the optical forces are non-conservative, the dynamical 
characteristics are very different from ordinary molecules 
bound by chemical forces; and the equilibrium configura- 
tions can correspond with either stable or a type of quasi- 
stable states in which the molecule maintains an average 
shape but individual particles exhibits periodic motion 
in the presence of frictional dissipation. Besides theo- 
retical interest, the concept of a photonic molecule may 
offer an alternative way to manipulate ultra-fine particles 
into artificial structures. We also note that light-induced 
forces in periodic photonic crystals has been calculated 
and discussed by Antonoyiannakis and Pendry 

Consider a cluster of N identical dielectric spheres 



illuminated by an incident time-harmonic EM field 
Einc{x,t) = Eqx exp{ikz — iut) + EqX exp(— ifcz — 
iujt), consisting of two counter-propagating x-polarized 
monochromatic beams of the same intensity, with angu- 
lar frequency uj and k = lo/c. The two beams interfere to 
form fringes with intensity varying as ~ cos(fcz), so that 
the gradient force along the z-direction traps the spheres 
in the xy-plane In order to be definitive, in what 
follows the dielectric spheres have radii rg = 0.414 fim, 
mass density 1050 kg/m'^, dielectric constant Er = 2.53 
(polystyrene), and the incident wavelength A = 0.52/xm 
{krs = 5). 

The intensity of the incident light is uniform in the 
plane, hence the forces on the spheres can only arise from 
multiple scattering. The light induced force (hence the 
strength of optical binding) is proportional to the inten- 
sity of the incident light, set to be 0.01 W//um^, similar 
to that used in the experiments in Ref. . We note that 
the physics of the problem remains qualitatively similar 
when the relevant parameters (sphere size, wavelength 
etc) are varied to within an order of magnitude, and a 
single plane wave is employed instead of two. In the latter 
all the phenomena are reproduced in the frame co-moving 
with the cluster along the z direction. 

We calculate the time-averaged electromagnetic force 
responsible for optical binding. The multiple scattering 
approach (MS) 9] is adopted to compute the electromag- 
netic fields, then the time-averaged force on the sphere 

i, < Fi >t is given by a surface integral of the time- 
averaged Maxwell stress tensor < T >t over the sphere's 
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MS-MST approach is highly accurate (retardation effect 
is included) and efficient, and it is the method of choice 
for a cluster of spheres in an arbitrary arrangement. 

Figure lA shows the force for a two-sphere cluster with 
the bi-sphere axis tilted at 45° to the polarization. The 
van der Waals force ^3 is included in the calculation. 
The alternating attractive-repulsive nature is a manifes- 
tation of the phase of the EM fields, and the radial force 
is seen to be stronger than the transverse force. The mag- 
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nitude of the force is inversely proportional to the spheres 
separation in the far zone, as the scattered radiation field 
decay inversely from the scatterer. This long range force, 
on the order of a pico-Newton, dominates over the short- 
range attractive van der Waals force at distances over a 
wavelength and oscillates between attraction and repul- 
sion with a period equal to the wavelength of the incident 
radiation. Stable radial positions are indicated by arrows 
in Fig. lA. They correspond to zero-force separations in 
which an increase in the radial separation would induce 
an attractive restoring force, whereas a reduction would 
induce a repulsive restoring force. 

The magnitude of the optical force depends on the ra- 
dius of the spheres. It can be fitted asymptotically to 
f{rs) cos[kR + (j){rs)]/R, where R is the separation be- 
tween the spheres, assumed to be S> Ts- For the po- 
larization perpendicular to the bi-sphere axis, is 
shown in Fig. IB. When fcr^ <C 1, the force is propor- 
tional to the sixth power of the radius. For kvg > 1 
the force exhibits complex variations, owing to resonant 
excitations. 

In Fig. lA, it is seen that the stable radial positions 
correspond to a positive aligning force along the trans- 
verse direction, therefore it is transversely unstable. The 
zero-force configurations can only be identified through a 
full two dimensional scan of the relevant forces, shown in 
Fig. IC. Multiple zero-force configurations exist on both 
the X and the y axes. For small separations, the force 
is much stronger when the bi-sphere axis is aligned with 
the polarization, due to the evanescent wave contribu- 
tion. For larger separations, the absence of propagating 
waves parallel to the polarization direction(s) means that 
the force is stronger when the bi-sphere axis is perpen- 
dicular to the polarization JJj . 

Zero force configurations are not necessarily stable. 
Perturbation calculations are required to locate stable 
equilibrium configurations. For general considerations, 
we assume a damping force proportional to the veloc- 
ity, as in most practical situations. Hence, for a clus- 
ter close to a zero force configuration denoted by the 

position vector — {xi,X2, ,X2i-i,X2i, ) 

|l3j, where {x2i-i,X2i) denotes the zero-force coordinate 
{x, y)i of the i*'' sphere, the linearized equation of mo- 
tion is m^^-^f- « KAx — b '^^^ , where m is the sphere 
mass, {/S.X2i-iT Ax2i) is the displacement vector of the 
i*^ sphere from the equilibrium, h the damping constant, 

and {K)jk = '^^'g^.Tk''' force constant matrix with 

flight being the light induced force. The eigenvalues of K 
dictate the stability of the cluster. Since a cluster bound 
by an incident light is an open system, the eigenvalues 
Xi can be conjugate pairs of complex numbers 14], as 
indeed found numerically with increasing probability as 
the number of spheres increases. 

Besides the two zero eigenvalues associated with 2D 




Separation / X 




Size Parameter q ^ |^ 

x/A 



FIG. 1: . (A) Forces on two dielectric spheres, plotted as a 
function of their separation. The bi-sphere axis is tilted at 
45'^to the polarization. The total radial force (optical plus 
van der Waals) and transverse force are shown by the black 
and red lines, respectively. Positive radial force means re- 
pulsion while a positive transverse force means the tendency 
to align the bi-sphere chain along the incident polarization. 
The green line denotes the radial force computed from the 
dipole approximation with radiative correction [Til . Signifi- 
cant difference in magnitude is seen. The thin blue line shows 
the van der Waals force. The arrows indicate the positions 
of stable radial separations. (B): [/(rs)] (see text) versus size 
parameter (krs). The spheres are aligned perpendicular to the 
polarization of the incident light and located at each other's 
far field zone. (C) Magnitude of the optical force on sphere 
2, with sphere 1 at the origin. The white region is forbidden 
(the spheres overlap). Blue, yellow and red corresponding to 
weak, medium and strong force respectively. Zero-force posi- 
tions are marked by the arrows, with green indicating stable 
positions and purple indicating unstable position. 



translation, if all the other eigenvalues are real and nega- 
tive, the cluster is stable. The appearance of real positive 
eigenvalue(s) or complex eigenvalues with Re{Xi) > de- 
notes instability. The eigenmodes for the complex eigen- 
values are associated with spiral motions. If the damp- 
ing coefficient b = 0, one of the modes spirals inwards 
upon perturbation, approaching the equilibrium position, 
while the other eigenmode spirals outward. If b is fi- 
nite, then Re{Xi) < indicates the existence of a critical 
damping constant bc„ucai = ^/rn\Im{Xi)\/ ^/\Re{Xi)\ at 
which the modes become stable for b > bcriUcai- We de- 
note this type of complex modes as quasi-stable modes. 
The appearance of one or more pairs of complex eigenval- 
ues with Re{Xi) < 0, with all the rest real and negative, 
denotes the cluster to be quasi-stable, i.e., potentially 
stable in the presence of damping. We note that the gra- 
dient trapping force in the z-direction is strong enough 
so that the system is stable if we consider the cluster as 
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a 3D object. 

For a multi-sphere cluster, the optical force can bind 
the spheres into a multiplicity of distinct geometries. In 
Fig. 2, the stable configurations are shown in panels (A)- 
(E), and quasi-stable configurations are shown in pan- 
els (F) and (G). For the particular stable configuration 
shown in Fig. 2A, the radial oscillation frequency is 540 
kHz (kT of energy, where T denotes room temperature, 
implies a fluctuation amplitude ~I4 nm), whereas the 
rotational oscillation frequency is 51 kHz {kT of energy 
implies a fluctuation amplitude ~0.09 radian). 

Comparison with molecules bound by chemical forces 
yields some interesting differences. First, atomic clusters 
can have multiple stable isomers, but the bond length al- 
ways has a narrow distribution. For photonic molecules, 
a fixed number of spheres can be bounded into a great 
number of stable configurations, each of different bound 
length and shape. Second, chemically bound molecules 
can rotate freely in space while photonic molecules, owing 
to the vector nature of light, are pinned in the two dimen- 
sional plane with a fixed orientation (relative to the light 
polarization). Third, for clusters without inversion sym- 
metry, photon scatterings are biased in one direction and 
hence the cluster absorbs light momentum, leading to 
"drifting" equilibrium states in which the cluster moves 
as a stable entity in the xy-plane. Thus configurations 
show in Fig. 2B and Fig. 2E are in drifting equilibrium 
while all the others are in static equilibrium. 




FIG. 2: . Examples of photonic molecules. Configurations 
(B) and (E) are in drifting equilibrium (see text), and the 
others are in static equilibrium. All eigenmodes of configu- 
rations (A)-(E) are stable, and (F)-(G) have either stable or 
quasi-stable modes. Size of the spheres is drawn to scale with 
respect to the wavelength. Polarization of the incident light 
beam is in the horizontal direction. (G) shows a 43-sphere 
photonic molecule. The red lines are drawn to unveil the 12- 
fold symmetric, square-triangle quasicrystal tiling. 

To study the dynamics of the quasi-stable photonic 
molecules, it is necessary to go beyond the linear stabil- 
ity analysis 0. We integrate the equations of motion 
using an adaptive time-step Runge-Kutta-Verner algo- 
rithm. The procedures are similar to molecular dynamics 
simulations, except that the inter-particle forces are now 
the light-induced forces. To illustrate the rich phenom- 
ena in dynamics, we study a six-sphere photonic molecule 
shown in Fig. 2F, which has all real negative eigenvalues 



except for one complex conjugate pair. We ignore van 
der Waals forces in the following calculations since the 
spheres are about 3 microns apart. At that distance the 
optical binding force is about 3000 times stronger than 
the van der Waals force. Figure 3 shows the cluster dy- 
namics at different levels of damping. For b > bcruicah 
each sphere exhibits damped oscillations that settle into 
the stable zero- force position (Fig 3D). For b slightly less 
than bcriticai, the zero- force position turns unstable and 
an attracting elliptic periodic orbit is formed (Fig 3C, pe- 
riod=15 fis) for each sphere in the immediate neighbour- 
hood of the zero force position. This periodic dynamics 
arise because the repulsive part of the force becomes less 
than that predicted by the linear stability analysis as the 
sphere spirals outward, allowing an orbit to be estab- 
lished. Here the damping dissipation is counterbalanced 
by the incident light energy. Hence from the dynam- 
ics perspective bcriticai is a supercritical Hopf bifurcation 
point [ig- A further decrease in b enlarges and deforms 
the periodic orbit (Fig 3B, period=32 /its), and additional 
decrease in damping leads to instability (Fig 3A). 

To investigate the thermal stability of the (stable) pho- 
tonic molecules, we move the spheres apart along the 
path of their lowest frequency eigenmode, until the pro- 
jection of the restoring force changes sign 16]. For typical 
photonic molecules the required dissociation energy per 
sphere was evaluated to be hundreds of kT {T= room 
temperature), e.g., the cluster in Fig. 2 A (two spheres 
at distance=3.16A) has a dissociation energy of 110 kT. 
Thus the photonic molecules are expected to be stable 
against thermal fluctuations. 

Although the input light field is homogeneous on the 
xy-plane, it is rather intriguing that the optical bind- 
ing force can arrange dielectric spheres into rather com- 
plex geometries — Fig. 2G shows a 43-sphere cluster 
with a quasicrystal-like geometry (red lines show the ex- 
act square-triangle quasicrystal tiling). The particular 
configuration shown has 18 pairs of complex eigenval- 
ues with their corresponding bcriticai between 1.3 and 
17 pN/ms-^ Thus if 5 > 17 pN/ms-\ a stable 
quasicrystal-like cluster can be realized. 

In summary, we showed using a rigorous multiple- 
scattering technique that the homogenous light-field of 
an incident plane wave can organize small particles into 
stable clusters. These clusters have well-defined geome- 
try (bond length and angle) and well-defined vibration 
frequencies, and they are bound by light, and hence we 
call them photonic molecules. Both the static and dy- 
namic properties are rather interesting, which is a man- 
ifestation that we have an open system, in contrast to 
chemical molecules are that bound by conservative forces. 
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FIG. 3: . Complex mode (centre-of-mass) trajectories at dif- 
ferent levels of damping for the cluster shown in Fig. 2F. (A): 
Trajectory of the spheres with 6 = 0, the cluster breaks up. 
(B)-(D) shows the trajectory of the selected sphere marked 
in (A). Black, red and blue colours respectively delineate the 
three consecutive time intervals of equal duration, from early 
to late. (B): For h = 3.71 pN/ms~^, in three 0.5 ms intervals. 
The cluster is attracted towards a periodic orbit. (C): For 
h = 4.31 pN/ms'^'bcriticai = 4.32 pN/ms"^), in three 130 ms 
intervals. The sphere is attracted towards an elliptic periodic 
orbit. (D): For h = 9.28 pN/ms"\ in three 0.2 ms intervals. 
The cluster is damped towards the stable zero-force position. 
The size of the periodic orbit is noted to be small compared 
to the interparticle separations. 
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